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Abstract 

We study the dynamics of hierarchy of piecewise maps generated by one-parameter 
famihes of trigonometric chaotic maps and one-parameter famihes of eUiptic chaotic 
maps of en and sn types, in detail. We calculate the Lyapunov exponent and Kolmogorov- 
Sinai entropy of the these maps with respect to control parameter. Non-ergodicity of 
these piecewise maps is proven analytically and investigated numerically . The invariant 
measure of these maps which are not equal to one or zero, appears to be characteristic 
of non-ergodicity behavior. A quantity of interest is the Kolmogorov- Sinai entropy, 
where for these maps are smaller than the sum of positive Lyapunov exponents and it 
confirms the non-ergodicity of the maps. 

Keywords: Non-ergodicity, piecewise maps, chaos, Lyapunov exponent, 

Kolmogorov- Sinai entropy. 
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Ergodicity theory is a branch of dynamic systems deaUng with questions of average. Some- 
times ergodic theory can make long-term predictions about the average behavior from initial 
data with hmited accuracy even for chaotic systems. The ergodic theory of chaos has been 
studied in detail by Eckmann and Ruelle [1, 2] and others, while non-ergodic mathematical 
models scarcely exist. It is shown that simplest one-dimensional dynamic systems satisfying 
the indecomposability assumption (and even the assumption of topological transitivity) may 
be non-ergodic, which shows that restrictions of this type are not quite reasonable in the 
context of general dynamical systems [3]. Also, ergodicity has not been proven for some 
systems and sometimes we have to face the problem of non-ergodicity of many chemical and 
physical systems especially in sohd systems (for example see [4, 5, 6, 7]). It is shown that, 
candidate solid state systems with extremely slow dynamics often have impurities remaining 
from their synthesis so there is always a suspicion that non-ergodicity comes from disorder 
[8]. Also, there are specific examples of arbitrary small perturbation to ergodic systems that 
are not ergodic [9, 10, 11]. The study of origins of the non-ergodicity and slow dynamics of 
polymer gels [12] and effects of temperature and swelling on chain dynamics in sol-gel tran- 
sition [13, 14] and also non-ergodicity transition in colloidal gelation [15] are some examples 
of research activity concerning non-ergodicity. 

The aim of present paper is twofold: To introduce the piecewise map with its invariant 
measure and to clarify its non-ergodic behavior. We present the hierarchy of one and many- 
parameter families of elliptic chaotic maps of cn and sn types which can generate non-ergodic 
behavior. In our definition of a piecewise map, we assume that it consist of components which 
can describe fixed point and chaotic behavior accordance with various values of parameter 
of the map. In particular, we argue that this map satisfy the non-ergodic assumption. The 
paper is organized as follows: The next section of the paper will note the competing defi- 
nitions of piecewise non-linear maps with complete boundary condition associated with the 
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one parameter families of chaotic maps and one-parameter families of elliptic chaotic maps 
of cn and sn types. In Section 3, we present the invariant measure of piecewise maps and 
in Section 4 the Kolmogorov- Sinai (KS) entropy and Lyapunov exponent (LE) of the map 
is studied. In Section 5, we review the ergodic theory and in Section 6, we give our results 
consist the KS-entropy and LE of the map and we shall explain how this piecewise maps 
behave non-ergodic. In addition there will be a concluding section and two appendices. 



2 Piecewise non-linear maps 

2.1 One-parameter families of trigonometric chaotic maps 

We first review the one-parameter chaotic maps which are used to construct the piecewise 
map. The one-parameter chaotic maps [16] are defined as the ratio of polynomials of degree 
N: 



di). . _ + 2Fi(-iV,iV,|,a;)) 



(a2 + l) + (a;2-l)(-l)iv ^F^{-N,N,'^,x) 



l + («2_i)(T^(^)2 



(2), , _ 2Fl(-iV,iV,i(l-^))) 



(«2 + 1) _ («2 _ l)(_l)iV ^p^^_N^ 1^ (1 _ ^)) 

2 



where N is an integer greater than one. Also, 

2F^{-N,N, ^,x) = (-l)^cos(2A^arccos V^) = (-l)^r2iv(v^) 

is the hypergeometric polynomials of degree N and Tiq{Un{x)) are Chebyshev polynomials 
of type I (type II), respectively. The conjugate maps of the one-parameter families of chaotic 
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maps which are used in derivation of their invariant measure and calculation of their KS- 
entropy are defined as: 



1 



(x, a) — ho (x,a) o h ^ — — tan^ (N arctan ^/x) 

1 1 

4>~^{x,a) — ho (f)^^\x,a) o h~^ — — cot^( A?" arctan —=) 

a"' ^yx 

Conjugacy means that invertible map h{x) = maps / = [0, 1] into [0, oo). To define the 
piecewise maps constructed from one-parameter chaotic map, we need to take into account 
boundary condition, namely we have to choose the states on the phase space. Here we 
present examples of these types which have been considered in the present paper: 



a\2x - 1)2 



Ax(l - x) + a^{2x - ly 



(2) _ Aa^xil - x) 

(P2 — 



1 + 4(q;2 - l)x(l - x) 



_ ^(2) 



a'^xi'ix - 3)2 



q;2x(4x-3)2 + (1-x)(4x-1)2 
Now, We introduce hierarchy of piecewise maps generated by one-parameter families of 
trigonometric chaotic maps. For piecewise map 0^ with even N, we have: 



(1). \ _ ) 



^^^[x,ai) aie[N,oo) 
{x, a2) 0L2 e [0, N\ 



The range of the parameters a\ and 0:2 in the maps ^^(a;, cvi) and ^^(a;, q;2) are cho- 
sen to grantee, respectively, chaotic behavior and two fixed points at x = and x — \. 
Figure 1 shows the plot of 02^''(x,a) for cti = 3 and «2 = 1- In ^^(•'J^! <^i) and ^^(a;, 0:2)! x 
is limited respectively to: x e [0.152, 0.848] and x e [0, 0.352] U [0.647, 1]. For given yo = 0.5 
the X are the roots of 0^ (x, a\) — yo and similarly the x are the roots of 0^ (x, 0:2) = yo- 
For piecewise map ^^(a;, a) with even N, we have: 
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The range of the parameters ai and a2 in the maps (p^ {x, cti) and (f)^ {x, a2) are chosen to 
grantee, respectively, two fixed points at x = and x — 1 and chaotic behavior. Figure 2 

(2) (2) (2^ 

shows the plot of 02 {x,a) for ai = 0.25 and a2 = 0.75. In (f)}^ {x,ai) and 0]^ (x, 0:2), x is 
limited respectively to: x e [0,0.378] U [0.621, 1] and x e [0.2,0.8]. For given yo = the x 

are the roots of (l)^^\x, cci) = yo and similarly the x are the roots of 0:2) = yo- 

(1 2) 

For piecewise map 0\r' with odd N, we have: 



This map has chaotic behavior for a2 and has a fixed point in a; = for ai. Figure 3 



X is limited respectively to: x e [0,0.2] U [0.315,1] and x e [0.086,0.453] U [0.947,1]. For 

(2) 

given yo = 0.5 the ij are the roots of (p)^ {x,ai) = yo and similarly the Xi are the roots of 
(I)'p{x,a2) = yo- 

2.2 One-parameter families of elliptic chaotic maps of cn and sn 



Here we first review a hierarchy of one-parameter families of elliptic of cn and sn types that 
have been used for constructing the piecewise maps with non-ergodic behavior. These kinds 
of maps are defined as the ratios of Jacobian elliptic functions of cn and sn types through 
the following equation [17]: 




a2 e [0, ^] U [TV, 00) 



shows the plot of (f)'^''^\x,a) for ai — 1.5 and a2 — 0.2. In (f)^^''^\x,ai) and (f)^^'^\x,a2), 



types 



0V(a;,a) 



a''(cn(Ncn-\^)))'' 



1 + (q;2 - l)(m(A^m-i(v^)))2 
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1 + (q;2 - l){sn{Nsn-^{y/x))y 



where a is control parameter. For N = 2 we have: 




q;2((1-A;2)(2x-1) + A;V)2 



(1 - F + 2k'^x - k^x^y + (q;2 - 1)((1 - k^){2x - 1) + k'^x^y 



(2), Aa'^x{l — k'^x){l — x) 



(1 - A;2x2)2 + Axil - x){a^ - 1)(1 - k'^x) 



It has been proved [17] that for small values of the parameter K of the elliptic function, 
these maps are topologically conjugate to the one parameter families of chaotic maps. Sim- 
ilar to the introduced piecewise in previous section (2.1), the piecewise of elliptic maps can 
be introduced. For an example for piecewise elliptic map 02 {^i «), we have: 



The range of the parameters ai and a2 in the maps 02 cti) and 02 (^2) are chosen to 
grantee, respectively, two fixed points at x = and x — 1 and chaotic behavior. Figure 4 
shows the plot of the elliptic map 02^^ (a;, a) for a\ — 0.5 and a2 — 2.5. In 02^^(a;, ai) and 
0f\a;,a2), x is limited respectively to: x E [0,0.28] U [0.72,1] and x G [0.027,0.973]. For 
given yo — 0.5 the x are the roots of 02^'* (x, cti) = yo and similarly the Xi are the roots of 
0f^(x,Q;2) = yo- 



Invariant measure or SRB measure is supported on an attractor and describes the statistical 
of long-time behavior of the orbits with respect to Lebesgue measure. For invariant measure 
of 0^'' map (i = 1,2) satisfying the Frobenius- Perron (FP) operator [18], we have: 




(l)y{x,ai) aiG[0, ^] 
(j)^2\x,a2) a2 e [^,00) 



3 Invariant measure 
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which is equivalent to: 

E /^(^)^ (3-1) 

For chaotic part of the piecewise map, i.e.; y e [0, yo] for 02^^(a;, ct) and y e [yo, 1] for both 
03^'^'*(a;, a) and 02^'* (x, a), the invariant measure iJ,{x,(3) is defined as: 

, ^ (3-2) 

With P > is the invariant measure of the maps (j)^^{x,a) provided that, we choose the 
parameter a in the foUowing form: 

r iV-l i 

El 2 I /^N iQ — k 

a = (3-3) 

in (x, a) maps for odd values of N, and; 

Q Y^l 2 J /^AT Q-k 

« = ^.E^T°^^^ (3-4) 

v^l 2 J /-lAT o-k 
2^k=0 ^2k+lP 

in (j)f^{x, a) maps for even values of N, where the symbol [ ] means greatest integer part. 
As it is shown in Appendix A, for satisfying the invariant measure ii{x, (5) with a, we obtain: 

with polynomials A{x) and B[x) defined as: 

[f] 

A{x) = C^,x^ (3-6) 



A;=0 



5(0:) = Cg+i^' (3-7) 

ifc=0 

If an invariant measure can be indecomposed into parts that are invariant, the measure is 
called non-ergodic. There may be several invariant measure for a dynamical system. If there 
is a fixed point x*, then a point distribution S{x — x*) in that point is an invariant measure. 
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even if the fixed point is unstable. Therefore for the fixed point part of the piecewise map 



average density measure ii{x, (5) has the following asymptotic form of the delta function as 
a goes to the zero and one, respectively; 



where the first one corresponds to invariant measure associated with the fixed point at x = 
and the latter one corresponds to fixed point Ski x — 1. 

Since for small values of K the parameter of the elliptic function, the elliptic chaotic maps of 
cn and sn types are topologically conjugate to the one-parameter families of trigonometric 
chaotic maps, we can obtain the invariant measure of these maps for small K [23]. As K 
vanishes, these maps are to trigonometric chaotic maps. 



4 Kolmogorov-Sinai Entropy and Lyapunov exponents 



Kolmogorov- Sinai (KS) entropy and Lyapunov characteristic exponents are two related ways 
of measuring 'disorder' in a dynamic system. The definition of them can be found in many 
textbooks [21]. To calculate KS-entropy here, we use the fact that it is equal to: 



which is also a statistical mechanical expression for the Lyapunov characteristic, that is 
mean divergence rate of two nearby orbits. As it is shown in Appendix B, the KS-entropy 
of ci;) following expression: 



i.e., y e [yo, 1] for 02^''(a;, a) and y e [0,yo] for both maps (t)'^''^\x,a) and 02^'* (x, a), the 







) 
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A useful numerical way to characterize chaotic phenomena in dynamic systems is by means of 
the Lyapunov exponents that describe the separation rate of systems whose initial conditions 
differ by a small perturbation. Suppose that there is a small change 5x{0) in the initial state 
x(0). At time t this has changed to Sx{t) given by: 

5x{t) ^ Sx{0)\^{x{0))\ = 6x{0)\<l)'{x{t - l)).0'(x(t - 2)) 0'(a;(O))|, (4-8) 

Where we have used the chain rule to expand the derivative of 4>. In the limit of infinitesi- 
mal perturbations dx{0) and infinite time we get an average exponential amplification, the 
Lyapunov exponent A , 

A = limt-,^j In \^^\ = Urrit^^j In |^(a;(0))| = lirrit-^^ ^ In (4-9) 

5 Ergodicity and non-ergodicity 

An probabilistic dynamic system is characterized as ergodic or non-ergodic by its marginal 
probability distributions. If the distributions have e.g. infinite variances so that a process 
mean cannot be defined, then the system is non-ergodic. An ergodic system has "convergent" 
qualities over time, variances are finite and a non-time-dependent process mean is clearly 
defined. Here a brief description of ergodic theory of chaos [19] is presented: let {fl, F, /i) 
be a probability space, Q, is the sample space, i.e.; the space of points, cu designating the 
elementary outcomes of an experiment. F is the o"-field (or a-algcbra) of events. An event 
is a set Acil which is of interest, the cr-field F is the ensemble of all events, i.e., AeF. 
Also fjL designates a probabihty measure of F. A transformation T is ergodic, if it has the 
probability that for almost every a;, the orbit {a;. To;, T^o;, ...} of a; is a sort of rephca of 
Q itself. Formally, we shall say that T is ergodic if each invariant set A, i.e.; a set such 
that T~^(A)=A, is trivial in the sense that it has measure either zero or one. T~^(A)=A 
=^//(A)=0 or ijl{A)—1. The transformation T is ergodic (or indecomposable, or metrically 
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transitive), if in the Birkhoff theorem, for any integrable, real- value function f, the limit 
value / is constant and we have /x-almost everywhere. 

In this case, the average value of /(.), evaluated along the orbit T^lu, converges /^-almost 
everywhere to the mathematical expecting or mean of /(.), evaluated on the space Q. In 
other word, for ergodic systems, the time average is equal to the space (or phase) average. 
One further consideration should be added at this point. The equality of KS-entropy and 
sum of all positive Lyapunov exponents; 

Ai>0 

indicates that in chaotic region, this map is ergodic as Birkhof ergodic theorem predicts [20]. 
In other words, when the KS-entropy is smaller than the sum of positive LE, the map has 
characterization of a non-ergodic behavior. 

Also, studied based on invariant measure analysis can be useful for confirming the non- 
ergodicity behavior of a map. For a non-ergodic system we have: 

li-\A) = {xe[0,l] y = M{x) ye A} 

fxiA) < 1 

i.e., the invariant measure which is not equal to zero or one, appears to be characteristic of 
non-ergodic behavior. 

6 Results and discussion 

In this section we present the results of the analysis numerically for piecewise map. Figure 
5, 6, 7 and 8 show the variation of LE and KS-entropy with the parameter a. A positive 
LE implies that two nearby trajectories exponentially diverge (at last locally). Negative LE 
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indicate contraction along certain directions, and zero LE indicate that along the relevant 
directions there is neither expansion nor contraction. 

In figure 5, the LE and the corresponding KS-entropy's have been shown for (l)^\x, ai) and 
02^^ (x, 02) by some points. A quantity of interest is that the KS-entropy is smaller than the 
sum of positive LE for piecewise maps. Because of this relation, it is clear that this map has 
characterization of a non-ergodic behavior. 

The invariant measure of this map is equal to 0.696 which is smaller than one, therefore this 
map behaves non-ergodic. The above analysis are presented for 02^"* (x, a), <j)^^''^\x^ a) and the 
elliptic map (l)^\x, a) (see figures 5-8 ).The invariant measure of these maps, respectively are 
equal to 0.6, 0.314, 0.946 which confirm the non-ergodicity behavior of introduced piecewise 
maps. 

7 Conclusion 

A recent attempts in introducing the hierarchy of chaotic maps with their invariant measure 
[16, 17, 22, 23, 24] allows us to advance in answering to a question how to define non-ergodic 
maps and what are the condition for non-ergodicity in these types of system. 
In this paper we introduce the piecewise maps with their invariant measure. Our numerically 
calculation shows that values of the KS-cntropy are smaller than the sum of positive LE for 
piecewise maps, therefore these maps behaves non-ergodic. Following, the non-ergodicity 
behavior, also can be confirmed by the invariant measure which is not equal to zero or one. 

8 Appendix A 

Similar to the calculation of the invariant measure in our pervious papers [16, 17, 22, 23, 24], 
we present here it for the piecewise chaotic map. In order to prove that measure (3-2) 
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satisfied equation (3-1), we consider the conjugate map; 



^l^Va;, a) = — tan^fiV arctan \/a;) (8"1) 



with measure /x^^ related to the measure n^^^ with the following relation: 

h^i^) = (TT^^'^^^iT^^- 

Denoting ^n{x, a) on the left-hand side of (8-1) by y and inverting it, we get: 

Xk = tan^(^ arctan ^yQ;2 + -^), k ^ 1, ...N. (8-2) 
Then, taking derivative of Xk with respect to y, we obtain 

Substituting the above result in equation (3-1), we get: 

+ "^y) = -^Yl + ^k)fi^^{xk)- (8-4) 

k 

Now, by considering the following ansatz for the invariant measure j^;p^{y)'- 
the above equation reduces to: 

1 + a^y _ o_ l+Xk . 
1 + Py ~ N^}l + PxJ 

which can be written as: 

1 + Q^y a (5 — 1 d 



T(in(nLi(r^+^fc)) )■ (8-6) 



To evaluate the second term in the right-hand side of above formulas, we can write the 
equation in the following form: 

= a^y cos^{N arctan y/x) — sm'^{N arctan y/x) 
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^ r'vi , rjv-i, . \ constant ^ 



where Xk are the roots of equation (8-1) and they are given by formula (8-2). Therefore, we 
have: 



dp 

d 



- In (1 - /3-i)^(Q;2y cos'^(N arctan yf^P^) - sin^N arctan yf^P^) ) 



NP , pN{l + a'y)Ail/P) 



P-1 iAil/P)yp^y + iBil/P)) 
In deriving the above formula, we have used the following identities: 

^, A(—x) . /,,- ^ Bi—x) 

cosfA' arctan Va^) = —-j^r, smfA^ arctan Va;) = yx-r^ 

(1 -I- xy^l'^ (1 -I- xy^i^ 

inserting the result (8-7) in (8-6), we get: 

1 + a^y 1 + a^y 



1 + Py B{l/P)/A{l/P) + P{aA{l/P)/B{l/P))y 
Hence, to get the final result, we have to choose the parameter a as: 

B{l/P) 



a 



9 Appendix B 



The KS entropy of one-parameter families of chaotic map is given by: 

d 



where 



/d 
li(x)dx In I —4>(x, a) I 
dx 



(p{x,a) = y = —{tan'^{N aTctan ^/x)) 



Therefore to calculate /i(//, (p{x, a)) we have: 

N 1 sin A'"( arctan ^/x) 



h{iJ,,ip{x,a)) = / iJ,{x)dx\n 







y/x{l + x) cos^ A'"(arctan ^/x) 
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using the relation (8-8), we get: 

N (l + x)^-^B(-x) 



. , 1 p ^/pdx 



(9-1) 



q;2 {A{-x)f 

We see that polynomials appearing in the numerator (denominator) of the integrand appear- 
ing on the right-hand side of equation (9-1), have ^[N — 1](|[A^]) simple roots, denoted by 
, k — 1, ^[N — l]{x^, k — 1, ^[N]) in the interval [0,oo). Hence we can writ the 
above formula in the following form: 

1 y^dx , /iv(i + x)^~inE-'^^'] 



Now making the following change of variabley^a^ = tan^, and taking into account that 
degree of numerators and denominator are equal for both even and odd values on Nwe get: 

1 /-oo jv 

h(n,^(x,a)) = - d0 {ln(— ) + (iV- l)ln|/? + l + (/3- l)cos^| 
TT Jo 

^^[{N-i)/2] In 1 1 _ + (1 + ^B^) COS e| - 3T.fJ? In 1 1 - + (1 + x^P) cos0\} 
Using the following integrals: 



1 

— In a -I- 6cos^ = < 

TT JO 



ln||| |a| < 1^1 
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Figures captions 

FIG. 1. plot of (f)2'\x, a) for ai — 3 and a2 — 1. 
FIG. 2. plot of (p^2\x, a) for ai = 0.25 and as = 0.75. 

(1 2) 

FIG. 3. plot of 03 ' {x, a) for ai = 1.5 and q;2 = 0.2. 

FIG. 4. plot of the elliptic map (f)^\x, a) for cci = 0.5 and a2 — 2.5. 

FIG. 5. The variation of the Lyapunov exponent (dotted curve) and the KS-entropy (solid 

curve) of the 02^'*(a;, a) with the parameter a. □ and -k show the valuates of Lyapunov ex- 
ponent and KS-entropy, respectively for (f)2\x,ai) and 02^'*(a;, q;2)- 

FIG. 6. The variation of the Lyapunov exponent (dotted curve) and the KS-entropy (solid 
curve) of the 02 {^^ ot) with the parameter a. □ and ★ show the valuates of Lyapunov ex- 

(2) (2) 

ponent and KS-entropy, respectively for 02 (^) cti) and 02 (^) (^2)- 

FIG. 7. The variation of the Lyapunov exponent (dotted curve) and the KS-entropy (solid 
curve) of the (f)^'^\x,a) with the parameter a. □ and ★ show the valuates of Lyapunov 
exponent and KS-entropy, respectively for (j)l^''^\x,ai) and 03^'^^ (x, 0:2). 
FIG. 8. The variation of the Lyapunov exponent (solid curve) of the elliptic map 02 {x, a) 
with the parameter a ior K — 0. □ and ★ show the valuates of Lyapunov exponent and 
KS-entropy, respectively for 02 {x, cti) and 02 {x, 0:2) • 
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